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Charge, Lattice, and Spin Dynamics in Photoinduced Phase Transitions from 
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To elucidate different photoinduced melting dynamics of charge orders observed in quasi-two- 
dimensional organic conductors fl-(BEDT-TTF) 2 RbZn(SCN) 4 and a-(BEDT-TTF) 2 I 3 [BEDT- 
TTF=bis(ethylenedithio)tetrathiafulvalene], we theoretically study photoinduced time evo- 
lution of charge and spin correlation functions on the basis of exact many-electron wave 
functions coupled with classical phonons in extended Peierls-Hubbard models on anisotropic 
triangular lattices. In both salts, the so-called horizontal-stripe charge order is stabilized 
by nearest-neighbor repulsive interactions and by electron-lattice interactions. In #-(BEDT- 
TTF)2RbZn(SCN)4 (abbreviated as #-RbZn), the stabilization energy due to lattice distortion 
is larger, so that larger quantity of energy needs to be absorbed for the melting of the charge 
and lattice orders. The photoinduced charge dynamics shows a complex behavior owing to a 
substantial number of nearly degenerate eigenstates involved. This is related to the high struc- 
tural symmetry when the lattice is undistorted. In Q-(BEDT-TTF)2l3 (abbreviated as a-h), 
the lattice stabilization energy is smaller, and smaller quantity of energy is sufficient to melt 
the charge and lattice orders leading to a metallic phase. The photoinduced charge dynamics 
shows a sinusoidal oscillation. In a-Is, the low structural symmetry ensures nearly spin-singlet 
bonds between hole-rich sites, where the spin correlation survives even after photoexcitation. 

KEYWORDS: photoinduced phase transition, charge order, metal-insulator transition, frustration 
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1. Introduction 

Electronic phases are controlled by changing pres- 
sure, constituents, or temperature in and near equilib- 
rium conditions. The effects of changing pressure or con- 
stituents can be introduced to theoretical calculations 
by changing model parameters. Electronic phases are 
dynamically controlled in nonequilibrium conditions as 
well, where the control process cannot be described by 
simply changing model parameters. Among them, pho- 
toinduced phase transitions have attracted much atten- 
tion. 1 ' 2 ) Along with the development of time-resolved 
spectroscopy, ultrafast transition dynamics enabled by 
strong electron correlation is one of the key issues at 
present. 

A considerable number of ultrafast photoinduced 
transitions toward metallic phases have been observed 
in molecular materials. Transitions from Mott-insulator 
phases arc realized in a halogen-bridged nickel-chain 
compound 3 ' 4 ^ and in a quasi-one-dimensional half- 
filled-band organic salt (BEDT-TTF) (F 2 TCNQ) 5 ' 
(F 2 TCNQ=2,5-dinuorotetracyanoquinodimethane). 
They are theoretically studied in refs. 6-9. Ultra- 
fast photoinduced transitions in charge density wave, 
Mott insulator, and metallic phases are investigated 
in an iodine-bridged platinum-chain compound. 10 ' 
Transitions from chargc-ordercd-insulator phases are 
realized in a quasi-onc-dimcnsional quarter-fillcd- 
band organic salt (EDO-TTF) 2 PF 6 11, 12 ' (EDO- 

*E-mail: kxy@ims.ac.jp 



TTF=cthylcncdioxy-tetrathiafulvalene) and in quasi- 
two-dimensional quartcr-fillcd-band organic salts 



0-(BEDT-TTF) 2 RbZn(SCN). 



13) 



-RbZn) and 



(BEDT-TTF) 2 I 3 13 ' 14) (a-I 3 ). The photoinduced dynam- 
ics in the former salt is also theoretically studied. 12, 15 ' 
The ground states of the latter salts are theoretically 
shown to be stabilized by nearest-neighbor repulsive 
interactions and by electron-lattice interactions. 16 ~ 20 ' 

The charge ordering in quasi-two-dimensional organic 
systems has been intensively studied both experimen- 
tally 21-23 -' and theoretically. 24, 25 ' The low-temperature 
phases are shown to have horizontal-stripe charge or- 
ders by NMR experiments (for 0-RbZn 26, 27 ' and for a- 
I 3 2 8,29)-j anc j thers. Their crystal structures are well 
analyzed by X-ray diffraction studies (for 6>-RbZn 30 ' 
and for a-I 3 31 '). Fi gure 1 schematically shows the low- 
temperature structures of the conduction layers with 
charge orders. The hole- rich sites indicated by the gray 
ellipses are aligned horizontally in this figure. For the 
optical properties of these salts, see ref. 22 for example. 
For (9-(BEDT-TTF) 2 CsZn(SCN) 4 (where Rb in 6>-RbZn 
is replaced by Cs) , which does not have a long-range or- 
der even in the nonmetallic phase at low temperatures, 
charge frustration effects are suggested to be important 
by permittivity measurements 32 -' and by the observation 
of thyristor behavior. 33 ' In #-RbZn, lattice effects are dis- 
cussed for the conductance and dielectric permittivity at 
various cooling and heating rates. 34 ' In CY-I3, ferroelec- 
tric polarization due to the charge ordering is recently 
reported, 35 ' where the ferroelectric order is of almost 
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Fig. 1. (Color online) Anisotropic triangular lattices for (a) 6- 
RbZn and (b) a-13. The rectangles denote the unit cells. The gray 
and white ellipses represent the hole-rich and hole-poor sites, 
respectively. The red or thickest lines possess the largest absolute 
values of transfer integrals after electron-lattice couplings are 
taken into account. 

purely electronic origin. 

After the pioneering theoretical work for the electronic 
states of quasi-two-dimensional organic conductors, 36 ' 
the charge-ordered states have been studied mainly in ex- 
tended Hubbard models because the nearest-neighbor re- 
pulsive interactions mainly stabilize them. 37-44 ) Charge 
frustration effects are recognized in various aspects. 45 ~ 48 ' 
Sometimes they are discussed in spinless fermion mod- 
els, 49,50 ' but the spin degrees of freedom are pointed out 
as important for determining the relative stability of dif- 
ferent charge ordering patterns. 19 ' Recent theories for 
charge ordcrings include refs. 51-55. In such frustrated 



systems, the lattice degrees of freedom are also impor- 
tant, 16 ~ 20 ' so that extended Peierls-Hubbard models on 
anisotropic triangular lattices are studied in the present 
work. 

As mentioned above, photoinduced melting of the 
charge orders has been observed in these salts by fem- 
tosecond reflection spectroscopy. 13 ' 14 ' It has revealed a 
fundamental difference between their dynamics, which 
would be due to different stabilities of these charge or- 
ders. The 6>-RbZn salt shows local melting and ultrafast 
recovery to the charge order in 0.2 ps irrespective of tem- 
perature and of excitation intensity. The CH-I3 salt shows 
critical slowing down in the transient reflectivity spec- 
tra, where decay times depend largely on temperature 
and on excitation intensity and vary from the picosec- 
ond timescale to the ten-nanosecond timescale. The rel- 
ative importance of nearest-neighbor repulsive interac- 
tions and electron-lattice interactions is shown to be dif- 
ferent between these two salts in our previous studies us- 
ing the unrestricted Hartree-Fock approximation, 18 ' the 
exact diagonalization, and the strong-coupling perturba- 
tion theory. 19 ' 

In #-RbZn, the lattice effect is strong and required to 
stabilize the horizontal-stripe charge order among differ- 
ent patterns of charge orders that are nearly degenerate 
if the lattice is undistorted and has the high structural 
symmetry. In a-I^, the lattice effect is weak because the 
lower structural symmetry (even without lattice distor- 
tion at high temperatures) already ensures some charge 
disproportionation. Namely, site B is hole-rich and site 
C is hole-poor in Fig. 1, which is regarded as a precursor 
to the horizontal-stripe charge order. Consequently, the 
electron-lattice coupling merely breaks the residual sym- 
metry at low temperatures. Therefore, different tenden- 
cies of these salts to the transition to a metallic phase are 
reasonably expected and preliminarily reported by using 
the time-dependent Hartree-Fock approximation. 56 ' 57 ' 

In the present study, we employ the exact many- 
electron wave functions coupled with classical phonons 
in order to elucidate the roles of electron-lattice interac- 
tions and different structural symmetries for the photoin- 
duced melting dynamics. Photoexcitation is introduced 
by an oscillating electric field, and the time evolution 
is obtained by solving the time-dependent Schrodingcr 
equation. 12,15 ' The absorbed energy that is required for 
the transition to a metallic phase is indeed larger for 0- 
RbZn than for C1-I3 indicating the higher stability of the 
charge order in 6*-RbZn. Different structural symmetries 
lead to different charge dynamics through different de- 
grees of degeneracy among charge ordering patterns. 

In reality, degrees of freedom outside the present model 
might affect the photoinduced dynamics. The differences 
between these salts are not limited to the conduction 
layers. Counter ions with different networks might in- 
fluence such nonlinear dynamics. In the present work, 
however, we focus on the conduction layers within two- 
dimensional models. Numerical results will be useful for 
understanding and for designing ultrafast and control- 
lable photoinduced dynamics in future. 
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2. Extended Peierls-Hubbard Model on Trian- 
gular Lattice 

For the quasi-two-dimensional organic conductors, we 
use the following two-dimensional 3/4- filled extended 
Hubbard model with electron-lattice couplings that mod- 
ulate transfer integrals on the anisotropic triangular lat- 
tice shown in Fig. i, 

H = ^ [(^±«iJ«ij)e i(e/ ' lc) * iJ - A(t) 4^ + H.c. 



t p 4 in f?-RbZn, 



and 



Ki 



L> 2uj 2 . ^ ' 



where c itT creates an electron with spin a at site 

i, nicr=c ia Cio ■> an d ni=J2 a n ia- The parameter ijj 
denotes the transfer integral for the bond between 
the neighboring ith and jth sites. Schematic illus- 
trations for the low-temperature structure of the 
conduction layer in 6>-RbZn 3 °) and that in a-13 31 ) 
are shown in Figs. 1(a) and 1(b), respectively. In 
this paper, we use the transfer integrals estimated 
from the extended Hiickel calculations based on 
the X-ray structural analyses 30, 31 - ) at high tempera- 
tures t itj =tfj as before, 17 ' 19 ^ tp , (=tg r =tg r )=0.035, 



HT 

P3 



-^ 4 T =0.095 for 



(9-RbZn and 



^=-0.0350, ^=-0.0461, ^ 3 T =0.0181, ^=0.1271, 
^ 2 T =0.1447, tgj T =0.0629, and ig^O.0245 for a-I 3 . 
Hereafter, we take eV as the unit of energy. The param- 
eter U denotes the on-site Coulomb repulsion strength, 
and Vi t j=V c or V p the intersite Coulomb repulsion 
strength between the neighboring ith and jth sites on 
the vertical or diagonal bonds, respectively. We apply the 
values [7=0.7, V c =0.35, and V p =0.3 to the interaction 
strengths for both salts, which reasonably reproduce the 
observed optical conductivity of f?-RbZn. 17 ) 

The quantity mj denotes the ith molecular displace- 
ment from the high-temperature structure, itjj its time 
derivative, ctij, Ki.j, and LUij the corresponding coupling 
constant, elastic coefficient, and bare phonon energy, re- 
spectively. The modulations of the transfer integrals from 
the high- to low-temperature phases are obtained by ap- 
propriate strengths of electron-lattice couplings, which 
reproduce the observed charge order with the help of 
Coulomb repulsion. 16 ~ 19 ) For the values of appropriate 
strengths, see ref. 17 for 6>-RbZn and ref. 19 for a-I 3 . 
Because each modulation is a monotonic function of the 
corresponding coupling strength, the latter is uniquely 
determined by setting the former at the experimentally 
observed value, i.e., the value estimated from the ex- 
tended Hiickel calculation based on the X-ray structural 
analysis at low temperature. The signs of (tjj ± ajjUjj) 
in the first line of eq. (1) are chosen so that yij > corre- 
sponds to the experimentally observed ones. For details, 
see ref. 17 for #-RbZn and ref. 19 for Q-I3. For simplicity 
in our present calculations, however, we take only the 
most important electron-lattice coupling, s^ 16 ) equiva- 
lent to s p 4 = ctp 4 / Kp^ 7 * 1 that makes t P 2 different from 



Ep2 



ip 4 T + u V AUp± 



tpi — tp4 CUp^Up^ , 

and si,2 = a^/ Kb2 18 ' 19 ' that makes t b 2 different from 
t b2 > in a-l 3 , 



tb2 — £fc2 T — a b2Ub2 



and 



ha 



if 2 + a- b2 u b2 



(1) 



Hereafter, the variables arc transformed as 



K, 



(2) 



for the bond index (i,j)=p4 in (9-RbZn and for (i,j)=b2 
in 0JT3. Now, we consider only one mode of molecular 
displacements. Then, we set the bare phonon energy at 



We use the 12-site clusters with periodic boundary 
conditions shown in Fig. 1. It should be noted that the 
12-sitc clusters arc tiled in a staggered manner to cover 
the anisotropic triangular lattice. As a consequence, the 
left-right symmetry is lost even if pl=p3 in 0-RbZn. The 
rectangle denotes a unit cell. The labeling of bonds cor- 
responds to the crystal structure at low temperatures. 
Hole-rich sites (A and B in 6>-RbZn and in a-13) are de- 
noted by gray ellipses, while hole-poor sites (C, D in 
#-RbZn and C, A' in a-l^) by white ellipses. Transfer 
integrals modulated by electron-lattice couplings are de- 
noted by different thickness. 

In the Pcicrls phase of the transfer integral, 
e i{e/hc)8 i<r A{t) ^ e denotes the absolute value of the elec- 
tronic charge, c the light velocity, dij the position vector 
Sij=Vj —Vi, and A(t) the time-dependent vector poten- 



tial, 



A(t) = -c [ At'E{t') , 
Jo 



(3) 



for the time-dependent electric field E(t). The field is 
here assumed to be given by 



E(t) = E ext 8(t)9(T irl . - t) sina; cxt i 



(4) 



with amplitude _E C xt, frequency u> G xt, and pulse width 
T , irr =27rA r cxt /tj cxt with integer iV ext . Here, 9{t) is the 
Heaviside step function, 8(t) = l for t >0 and 6*(t)=0 for 
t <0. By setting the intermolecular distance along the 
c-axis of 0-RbZn and that along the a-axis of 0J-I3 as the 
unit of length and e at unity, the amplitude of the electric 
field I .Eoxt I is given in units of eV. 

The time-dependent Schrodinger equation for the ex- 
act many-electron wave function | ip(t)) is numerically 
solved by 



ip(t + dt)) ~ exp 



-dm 

h 



dt 



\m), (5) 
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Fig. 2. (Color online) Optical conductivity spectra obtained 
by exact diagonalization of the 12-site clusters with t/=0.7, 
V c =0.35, and V p =0.3 for (a) 0-RbZn (s p4 =0.1) and (b) «-I 3 
(s;,2=0.1) with different polarizations as indicated. The peak- 
broadening parameter rj is set at 0.005. 



where the time evolution operator is expanded as 



cxp 



-dtH t + 



dt 



-dtH t + 



dt 



(6) 

with time slice cft=0.054 for 6»-RbZn and di=0.050 for a- 
I3 up to the order by which the deviation of (ip(t) \ 
from unity is less than 10~ 12 . The initial state \ip(0)) is 
set at the ground state. The classical equation of phonon 
motion is solved by the leapfrog method, where the force 
is derived from the Hellmann-Fcynman theorem. The ini- 
tial displacement y(0) is at the minimum of the adiabatic 
potential for the ground state, so that the initial velocity 
?/(0) is zero. 

3. Conductivity Spectra 

Before proceeding to the photoinduccd dynamics, we 
calculate the regular parts of the optical conductivity 
spectra by the formula, 



a(uj) = 



1 

No;' 



-Im 



(4 



1 3- 



1 



En — H 



j I V>0> 



(7) 

and show them in Fig. 2. Here, IV'o) denotes the ground 
state, E its energy, j the current operator with rigid 
transfer integrals, 77 the peak-broadening parameter, and 



N the system size. Namely, the infrared conductivity 
originating from phonons is ignored. These spectra corre- 
sponding to w p h=0 are useful for judging whether or not 
the frequency of the oscillating electric field (the energy 
of the pump light) is resonant with any charge-transfer 
excitation. In both cases of 6>-RbZn and of a-l3, the oscil- 
lator strengths are larger for polarization parallel to the 
stripes (a- and 6-axes, respectively) than for polariza- 
tion perpendicular to them (c- and a-axes, respectively). 
This is because the transfer integrals are larger for charge 
transfers along the p- and 6-bonds, respectively, than 
those along the c- and a-bonds, respectively 

The charge-transfer excitation with the largest oscilla- 
tor strength in ^-RbZn [Fig. 2(a)] has a lower energy and 
a smaller oscillator strength than that in a-13 [Fig. 2(b)]. 
For comparison, we will use similar excitation energies 
w ex t for the two compounds, which are almost resonant 
with excited states with substantial oscillator strengths. 
The targeted excited state with w=0.2903 in 0-RbZn has 
charge densities, 2 — (n^), 0.67 at sites A and B, and 
0.33 at sites C and D, while that with cl>=0.3162 in a-^ 
has 0.693 at site A, 0.648 at site B, 0.326 at site C, and 
0.333 at site A'. Note that the ground state in 6*-RbZn 
has charge densities, 0.93 at sites A and B, and 0.07 at 
sites C and D, while that in a-^ has 0.91 at site A, 0.78 at 
site B, 0.14 at site C, and 0.17 at site A'. Because of the 
smallness of the present clusters, these values themselves 
are insignificant, but they are useful when we analyze 
the respective photoinduced dynamics later. Concerning 
the overall shapes of the conductivity spectra, the time- 
dependent Hartree-Fock approximation for much larger 
system sizes gives smaller differences between <j(lu) for 
6*-RbZn and ct(lu) for a-lz than shown here. 58 ' The ex- 
perimentally observed spectra show very broad features 
for both salts, 22 ' so that a comparison is not straightfor- 
ward. 

4. Photoinduced Melting Dynamics 

For photoexcitations hereafter, we use the polarization 
parallel to the stripes and similar energies that are nearly 
resonant with a charge-transfer excitation with substan- 
tial oscillator strength. The action of photoexcitations is 
basically to destroy the charge order, so that excitation- 
dcnsity-dcpcndcnt quantities arc insensitive to the po- 
larization and to the frequency of the oscillating electric 
field. 15 ' 

The time evolution of the charge densities, 2 — (rii), 
during (t < Ti rr ) and after (t > T- llr ) photoexcitation is 
shown in Fig. 3. Because we take eV as the unit of energy, 
t=1520 corresponds to 1 ps. In both cases of 0-RbZn and 
of CH-I3, the amplitudes of the oscillating electric fields 
are chosen so that they are near critical values for melt- 
ing the charge orders at t=T- m with w ox t shown in the 
caption and with iV cx t=10. The total energy per site is 
increased by 0.03152 for 6>-RbZn and by 0.00963 for a- 
I3. The resonantly-excited charge-transfer processes are 
clearly visible, which have periods that are nearly equal 
to those of the oscillating fields 27r/o; ex t [21.6 in Fig. 3(a) 
and 19.9 in Fig. 3(b)]. In the figure, the time dependence 
of the electric field, eq. (4), is shown by the thin lines for 
comparison. The amplitudes of such oscillating charge 
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Fig. 3. (Color online) Time dependence of charge densities dur- 
ing (t < Ti rr ) and after (t > Tj rr ) photoexcitation along the 
stripes for (a) 0-RbZn (w oxt =0.2905, £ cx t=0.048) and (b) q-I 3 
(aj ex t=0.3165, i? C xt =0.007). The parameters shown in Fig. 2 and 
w ph = 0.01 are used. The thin lines show the time dependence of 
the electric field, eq. (4). Its period, 27r/cj ex t, is 21.6 in (a) and 
19.9 in (b). 

densities with these frequencies are slightly larger for a- 
I3 . There is a nonnegligible difference in the overall time 
dependences. The oscillation for 6*-RbZn is more complex 
[Fig. 3(a)], while that for o>±3 consists of a fewer com- 
ponents [Fig. 3(b)]. We always find such difference irre- 
spective of polarization or energy of the photoexcitation 
(not shown). Such difference is also seen by the time- 
dependent Hartrec-Fock approximation for much larger 
system sizes, 58 ' so that systematic differences seem to ex- 
ist in the electronic states. The origins of these behaviors 
are discussed below. 

To analyze the character of the electronic wave func- 
tion at time t, \ tp{t)), we can expand it by using some 
orthonormal functions. Because the applied field is uni- 
form, the molecular displacements yij (t) are homoge- 
neous, Ui,j(t) = y{t), which is, of course, a classical 
variable and is obtained as the solution to the classi- 
cal equation of phonon motion. For each y(t), we em- 
ploy the cigenfunctions of the Hamiltonian, | ipjfy{t)]), 
for such orthonormal functions. At t—Q, \ ip(0)) is set at 
the ground state | ipo[y(0)}). At t >0, | ip(t)) deviates 
from the lowest-energy state | ipo[y(t)]) of the Hamilto- 
nian with y(t). Then, let us define the ground-state pop- 



Fig. 4. (Color online) Time dependence of ground-state popula- 
tions I (ipo[y{t)] I 1 2 an d modulations of transfer integrals y 
(inset) during and after photoexcitation along the stripes for (a) 
6-RbZn with y=y p 4 and (b) a-Iz with y=yb2- The parameters 
shown in Fig. 3 are used. 

ulation by the square of the absolute value of the over- 
lap between the wave function, | ip(t)), and the lowest- 
energy state I ipo[y(t)]) for the molecular displacement 
y(t). Technically, we obtain this quantity by projecting 
the wave function | ip(t)) to | ipo[y(t)]) at each step of 
time evolution. If the time evolution is adiabatic, this 
quantity is independent of time. 

The time dependence of the ground-state population 
is shown in Fig. 4. The pulse width T m is denoted by 
the triangles. In 0-RbZn, the ground-state population 
is about 0.4 at i=Xi rr and shows a complex behavior 
for t > Tj rr [Fig. 4(a)]. It varies significantly when the 
displacement y p i{t) is small. As discussed in the follow- 
ing paragraph, the undistortcd {yij—0) lattice structure 
possesses the high symmetry in 0-RbZn, leading to quasi- 
degeneracy of excited states with different charge distri- 
butions. The weights of these excited states in \ tp{t)) are 
sensitively altered by y p i(t). This causes the variation of 
the ground-state population near y p ±(t)=0. In a-13, on 
the other hand, the ground-state population is about 0.65 
at i=Xi rr , which is larger than in 0-RbZn, and becomes 
independent of time for t > T lrl [Fig. 4(b)]. In contrast 
to 0-RbZn, the undistorted (?/ij =0) lattice structure pos- 
sesses quite low symmetry in a-l3 , which is free from such 
(quasi-)dcgeneracy. This causes a constant ground-state 
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Fig. 5. (Color online) Adiabatic potentials of ground and excited 
states for (a) 0-RbZn and (b) a-13. Those of the resonantly- 
excited charge-transfer states are presented in blue (by thick 
lines). The trajectories of the electronic energies with y p 4(t) and 
Vb2 ft) are plotted in red (by crosses) for the parameters used in 
Fig. 3. 

population for t > T; rr . The ground-state populations at 
t=T- lxt are rather insensitive to the polarization and to 
the energy of the photoexcitation. For any combination 
of the polarization and the energy, the ground-state pop- 
ulation at t=T[ TT for 0-RbZn is smaller than that for CK-I3. 
This indicates that a smaller number of excited states are 
involved in the dynamics for a-13. The oscillations of the 
charge densities shown in Fig. 3(b) are nearly sinusoidal 
because only a few excited states are involved in the pho- 
toinduccd dynamics for a-13. Both of the displacements 
y p i and yvi show an oscillator's behavior in the insets of 
Fig. 4. The periods of the oscillations are about 1500 to 
2000, which are much larger than those of bare phonons 
(27r/w p h=628) owing to the electron-lattice interaction 
(phonon softening). The oscillators seem to be weakly 
damped, but they do not vanish within accessible com- 
putation times. 

For the sake of discussions later, let us set aside the 
time evolution for the moment and calculate the total 
energies for the eigenstates | ipj [y(t)]) on the assump- 
tion of y = 0. The adiabatic potentials for the ground 
and excited states obtained in this manner are shown in 
Fig. 5, as a function of the molecular displacement y P 4 
or yb2- For both 0-RbZn and a-13, there are resonantly- 



excited charge-transfer states in Figs. 3 and 4, which 
have substantial oscillator strengths and excitation en- 
ergies of about uj cx t ~0.3 (Fig. 2) at the potential min- 
ima for the ground states (i.e., at y c±0.05). They are 
indicated in blue (by thick lines). For excited states 
with much higher energies than those of the resonantly- 
excited charge-transfer states, their adiabatic potentials 
are omitted for the sake of clarity. First of all, the lattice 
stabilization energy is larger for 6*-RbZn [0.114 in total or 
0.0095 per site in Fig. 5(a)] than for a-I 3 [0.019 in total 
or 0.0016 per site in Fig. 5(b)]. In 0-RbZn, the energies 
are distributed more densely at y P 4=0 than in a-13. It is 
because different patterns of charge orders are nearly de- 
generate owing to the high symmetry of its undistortcd 
lattice structure and the resultant charge frustration. 17 ^ 

Now, we return to the analysis of the photoinduced 
dynamics. Here we show the trajectory of the electronic 
energy, which is obtained by subtracting the lattice ki- 
netic energy, Yliij) (Kij/^ij^ij' ^ rom t ^ ie tota ^ ener §y- 
It should be noted again that, for both 6>-RbZn and a-^, 
the amplitudes of the oscillating electric fields, E cxt , are 
chosen so that they are near critical values for melting 
the charge orders. The total energy is of course conserved 
after photoexcitation (t > T u . r ). The electronic energy 
and the lattice kinetic energy are not monotonic func- 
tions of time because some energy is transferred back 
and forth. The electronic energy after photoexcitation 
(t > Ti rr ) is not a symmetric function of y. It is clearly 
seen for 0-RbZn in Fig. 5(a), where the ground-state pop- 
ulation shows a complex behavior [Fig. 4(a)]. Quantita- 
tively, the lattice kinetic energy is small in the present 
model, which restricts the number of phonon modes to 
one. In general, with an increasing number of phonon 
modes that are coupled to the electronic system, more 
energy would be transferred to the phonon system. How- 
ever, the electronic energy dominates the total energy 
because the nearest-neighbor repulsive interactions are 
much larger than the electron-lattice interactions. 

The electronic energy reaches an energy much higher 
than that of the resonantly-excited charge-transfer state 
for (9-RbZn [Fig. 5(a)]. This indicates that quite a large 
number of excited states are involved in the photoin- 
duced dynamics. In a-13, on the other hand, the energies 
are distributed sparsely around j/b2=0 owing to the low 
symmetry even in the undistortcd structure. Recall that 
charge disproportionation already exists between sites B 
and C at j/b2=0, and a finite y^ breaks the residual sym- 
metry between sites A and A'. The trajectory of the elec- 
tronic energy reaches an energy that is much lower than 
that of the resonantly-excited charge-transfer state for 
a-^ [Fig. 5(b)]. This indicates that a much smaller num- 
ber of excited states are involved in this photoinduced 
dynamics. 

5. Excitation-Density-Dependent Correlations 

Hereafter, charge densities and correlation functions 
arc calculated with different amplitudes, E cxt - For each 
E'ext, they are averaged after t=T- lrv over the period of 
r irr =27riVe X t/wext < t < 2irN ohs /Lu cxt with A cxt =10 and 
Aobs=100. They are shown as a function of the increment 
in the total energy per site AE. 
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Fig. 6. (Color online) Time averages of charge densities and 
modulations of transfer integrals (inset) after photoexcitation 
along the stripes for (a) 6-RbZn (a; ex t=0.2905) and (b) ct-13 
(aj ex t =0.3165), as a function of increment in total energy AE. 
The parameters shown in Fig. 2 and oj p i 1 =0.01 are used. 
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Fig. 7. (Color online) Time averages of double occupancies after 
photoexcitation along the stripes for (a) 0-RbZn (u e xt=0.2905) 
and (b) a-13 (u e xt =0.3165), as a function of increment in total 
energy AE. AE C is indicated by the arrows. The parameters 
shown in Fig. 2 and oj p i 1 =0.01 are used. 



The charge densities, 2 — (rij) = (nj 1 ) with n^=J2 a n^, 

n i<j= c ia c icn an d c 'i a — C L ' at f° ur sites m the unit cell of 
Fig. 1 after photoexcitation are plotted in Fig. 6. In each 
case, there is a critical value of the increment in the total 
energy per site AE Cl above which the interaction- induced 
charge disproportionation disappears. The critical value 
[A£; c =0.031 in Fig. 6(a)] for 6>-RbZn is larger than that 
[AE^O.OIO in Fig. 6(b)] for a-I 3 . This is understandable 
in light of the fact that the lattice stabilization energy is 
larger for 0-RbZn. We have calculated the critical values 
with different polarizations and energies of photoexcita- 
tions. For the polarization perpendicular to the stripes 
and for the same excitation energies as in this section, 
the critical value for #-RbZn is A£' c =0.020 and that for 
CH-I3 is A£? c =0.008. For the polarization parallel to the 
stripes and for cj oxt =0.25, which is far from any charge- 
transfer excitation energy (Fig. 2), the critical value for 
0-RbZn is A£ c =0.028 and that for a-I 3 is A£ c =0.014. 
We always find the relation, AE C for 6>-RbZn is larger 
than AE C for a-13. 

Above AE C , the time- averaged charge densities are 
uniform for f?-RbZn, but this is not the case for a-l 3 with 
low structural symmetry. For a-13, the time-averaged 
charge density at site A and that at site A' are equal 
above AE C . The difference between those at sites B and C 



decreases with increasing photoabsorption, and it finally 
vanishes after photoabsorption much larger than AE C . 
The insets of Fig. 6 show that the molecular displace- 
ments that modulate the transfer integrals also vanish 
around AE C . In both cases, the electronic states above 
AE C are regarded as metallic because of almost vanishing 
charge correlations shown later. 

The double occupancies, (n^nh), after photoexcita- 
tion are plotted in Fig. 7. Their behavior is quite similar 
to that of the charge densities. Above the same criti- 
cal values of AE C , the interaction- induced differences in 
the double occupancies disappear. In 0-RbZn, the dou- 
ble occupancy is large at hole-rich sites and very small 
at hole-poor sites below AE C [Fig. 7(a)]. The difference 
becomes small with increasing photoabsorption and dis- 
appears at AE C . Above AE C , the time- averaged double 
occupancies are uniform and increase with photoabsorp- 
tion for £?-RbZn. For a-13, that at site A and that at site 
A' are equal above AE C [Fig. 7(b)]. The difference be- 
tween those at sites B and C decreases with increasing 
photoabsorption, and it finally vanishes after photoab- 
sorption much larger than AE C . In both cases, the dou- 
ble occupancies approach 0.0625=(l/4) 2 for very large 
photoabsorption, which is the value the noninteracting 
uniform state possesses. 
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Fig. 8. (Color online) Time averages of charge-charge correlation 
functions between hole-rich and hole-poor sites after photoex- 
citation along the stripes, (a-b) for S-RbZn (u ext =0.2905) and 
(c-d) for a-13 (u) e xt =0.3165), as a function of increment in to- 
tal energy AE. AB C is indicated by the arrows. The parameters 
shown in Fig. 2 and tj ph =0.01 are used. 



Fig. 9. (Color online) Time averages of (a,c) charge-charge cor- 
relation functions and (b,d) spin-spin correlation functions be- 
tween hole-rich sites after photoexcitation along the stripes, (a-b) 
for 0-RbZn (u%xt =0.2905) and (c-d) for a-l 3 (w cxt =0.3165), as 
a function of increment in total energy AE. AE C is indicated by 
the arrows. The parameters shown in Fig. 2 and oj ph =0.01 are 
used. 



The charge-charge correlation functions, (n^n^), with 
i=l on (hole-rich) site A and j=4, 5, or 6 on (hole-poor) 
site C [Figs. 8(a) and 8(c)] or j=7, 8, or 9 on (hole-poor) 
site D [Fig. 8(b)] or A' [Fig. 8(d)] after photoexcitation 
are plotted here for 0-RbZn [Figs. 8(a) and 8(b)] and for 
CV-I3 [Figs. 8(c) and 8(d)]. Those between the nearest- 
neighbor sites are indicated in the graph legends by the 
labeling of bonds. These charge-charge correlation func- 
tions between hole-rich and hole-poor sites are small be- 
fore photoexcitation (they are zero in the strong-coupling 
limit). They increase with photoabsorption below AE C , 
and change only a little above AE C . Between the dis- 
tant pair of sites [for j=5 in Figs. 8(a) and 8(c), and for 
j=9 in Figs. 8(b) and 8(d)], it is near 0.25=(l/2) 2 above 
A.E C , which is the value the noninteracting uniform state 
possesses. In other words, the charge-charge correlations 
almost disappear above AE C , indicating that these states 
are metallic with photocarriers. 

The charge-charge correlation functions with i=l 
and .7=10, 11, or 12 between hole-rich sites A and B 
[Figs. 9(a) and 9(c)] and the spin-spin correlation func- 
tions, {S*Sj) with Sf — (rii| — n%i)/2, for the same 
combinations of i=l and j [Figs. 9(b) and 9(d)] after 
photoexcitation are plotted here for 6>-RbZn [Figs. 9(a) 
and 9(b)] and for a-13 [Figs. 9(c) and 9(d)]. The corre- 
lation functions between the nearest-neighbor sites are 



indicated in the graph legends by the labeling of bonds. 
The charge-charge correlation functions between hole- 
rich sites are large before photoexcitation (they are unity 
in the strong-coupling limit). They decrease with increas- 
ing photoabsorption below AE C , and change only a little 
above AE C . Between the distant pair of sites and 
J=12), it is near 0.25=(l/2) 2 above AE C [Figs. 9(a) and 
9(c)], which is the value the noninteracting uniform state 
possesses. 

For 0-RbZn, the spin-spin correlation functions are 
negative and large on the p4 bonds (the difference be- 
tween j=10 and j=ll is due to the loss of the left-right 
symmetry in the present tiling of the 12-site clusters) 
before photoexcitation [Fig. 9(b)]. Their magnitudes de- 
crease with increasing photoabsorption below AE C , and 
they are very small and almost independent of the pair of 
sites above AE C . Between the distant pair of sites (i=l 
and j=12), it is always very small. For a-La, the spin- 
spin correlation functions are significant only on the 62' 
bonds before photoexcitation [Fig. 9(d)]. Its magnitude 
decreases with increasing photoabsorption below AE C , 
but it remains finite just above AE C due to the large 
tb2 even without lattice distortion (though it becomes 
very small for AE much larger than AE C ). The other 
spin-spin correlation functions are always very small not 
only between the distant pair of sites but also on the 63 
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bonds. Just above AE C , the spin-spin correlations almost 
disappear for 0-RbZn [Fig. 9(b)], but they remain nearly 
spin-singlets on the b2' bonds for a-±3 [Fig. 9(d)]. This 
can be interpreted as follows. The energy transfer from 
the oscillating electric field to the spin degrees of free- 
dom in the electronic state does (does not) reach 100% 
in 0-RbZn (a-I 3 ). 

6. Summary 

Differences between the photoinduccd melting dynam- 
ics of the charge order in 6>-(BEDT-TTF) 2 RbZn(SCN) 4 
(0-RbZn) and that in a-(BEDT-TTF) 2 I 3 (a-I 3 ) are 
studied in two-dimensional 3/4- filled extended Peierls- 
Hubbard models on anisotropic triangular lattices. The 
charge ordering patterns are very similar and known 
as horizontal stripes. They are stabilized by nearest- 
neighbor repulsive interactions and by electron-lattice 
interactions, but their relative importance in stabilizing 
the charge orders is different between these salts. 16 ~ 20 ) 
Time evolutions of exact many-electron wave functions 
coupled with classical phonons are obtained by solving 
the time-dependent Schrodingcr equation during and af- 
ter the application of an oscillating electric field. Cal- 
culations are performed for photocxcitations with differ- 
ent polarizations and energies. The qualitative difference 
between these dynamics is basically consistent with the 
experimentally observed one: 13 ) the photoabsorption re- 
quired for the transition to a metallic phase is larger for 
0-RbZn than for a-I^ indicating the higher stability of 
the charge order in 0-RbZn. 

Although the charge ordering patterns are very simi- 
lar, these salts have different crystal structures. The high 
symmetry in the high-temperature structure of (9-RbZn 
causes charge frustration leading to quasi-degeneracy 
among electronic states possessing different charge or- 
dering patterns if lattice distortions are absent. When 
the lattice distortions are small, the energy distribution 
is so dense that a large number of excited states are 
involved in the photoinduced dynamics. This causes a 
complex time dependence of the charge densities. The 
low symmetry in the high-temperature structure of a- 
I3 causes a precursor to the charge order without lat- 
tice distortion. The energy distribution is rather sparse, 
and a fewer number of excited states are involved in the 
photoinduced dynamics. This causes a sinusoidal time 
dependence of the charge densities. 

Charge-charge and spin-spin correlation functions as 
well as charge densities and molecular displacements are 
averaged over time after photoexcitation and plotted as 
a function of the increment in the total energy per site 
AE. There is a critical value AE C in the increment, above 
which the interaction-induced charge disproportionation 
and lattice distortion disappear. The critical value AE C 
is larger for 0-RbZn. This relation always holds irrespec- 
tive of polarization or energy of the photoexcitation. The 
correlation functions monotonically change with pho- 
toabsorption toward AE C from below. Their changes are 
small above AE Cl where the charge-charge correlations 
almost disappear and the system is regarded as metallic. 
Only on the bonds between hole-rich sites, the spin-spin 
correlations can be large before photoexcitation. Above 



AE Cl they become very weak and are still non-negligible 
only on the 62' bonds in CK-I3, where the transfer inte- 
grals are large even without lattice distortion owing to 
the low structural symmetry. 

From the theoretical viewpoint, the difference in the 
photoinduced dynamics of these salts is a quantitative 
one. The higher stability of the charge order in 0-RbZn is 
due to large stabilization energy. In reality, lattice anhar- 
monicity might be important. For instance, the difference 
in the insulating layers of counter ions [RbZn(SCN)4 in 
0-RbZn and I3 in a-13] might also contribute to the sta- 
bility of the charge order in (9-RbZn compared with that 
in CK-I3. What is peculiar to a-^ is such low structural 
symmetry that charge disproportionation exists between 
sites B and C even in the metallic phase and the transi- 
tion to the metallic phase is achieved by dissolving only 
the charge disproportionation between sites A and A'. 
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